Competition between Spin-Orbit Interaction and Zeeman Coupling in Rashba 2DEGs 
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We investigate systematically how the interplay between Rashba spin-orbit interaction and Zee- 
man coupling affects the electron transport and the spin dynamics in InGaAs-based 2D electron 
gases. From the quantitative analysis of the magnetoconductance, measured in the presence of an 
in-plane magnetic field, we conclude that this interplay results in a spin-induced breaking of time 
reversal symmetry and in an enhancement of the spin relaxation time. Both effects, due to a par- 
tial alignment of the electron spin along the applied magnetic field, are found to be in excellent 
agreement with recent theoretical predictions. 

PACS numbers: 73.23.-b, 71.70.Ej, 72.25.Rb 



Achieving control of the orbital motion of electrons by 
acting on their spin is a key concept in modern spin- 
tronics and is at the basis of many proposals in the field 
of quantum information[l|. Two physical mechanisms 
are used to influence the dynamics of the electron spin 
in normal conductors: spin-orbit interaction (SOI) and 
Zeeman coupling. In the presence of elastic scattering, 
these two mechanisms affect the spin in different ways. 
SOI is responsible for the randomization of the spin direc- 
tion whereas the Zeeman coupling tends to align the spin 
along the applied magnetic field. Depending on the rela- 
tive strength of these interactions, this interplay of SOI 
and Zeeman coupling is responsible for the occurrence of 
a variety of physical phenomenal^ |3| . 

Quantum wells (QW) that define 2-dimensional elec- 
tron gases (2DEGs) are particularly suitable for the ex- 
perimental investigation of the competition between SOI 
and Zeeman coupling, since they give control over many 
of the relevant physical parameters. Specifically, in these 
systems the SOI strength can be controlled by an appro- 
priate QW designjjl and by applying a voltage to a gate 
electrode |5|- The electron mobility is usually density de- 
pendent, so that the elastic scattering time can also be 
tuned by acting on the gate. Finally, Zeeman coupling 
to the spin can be achieved with minimal coupling to the 
orbital motion of the electrons by applying a magnetic 
field parallel to the conduction plane. 

In this Letter we study the competition of SOI and 
Zeeman coupling via magnetoconductance measurements 
in InGaAs-based 2DEGs with different Rashba SOI 
strength. From the detailed quantitative analysis of the 
weak antilocalization as a function of an applied in-plane 
magnetic field (-Bii), we find that the partial alignment of 
the spin along B\\ results in a spin-induced time reversal 
symmetry (TRS) breaking, and in an increase of the spin 
relaxation time. The increase in spin relaxation time is 
found to be quadratic with B\\, and strongly dependent 
on the SOI strength and the elastic scattering time. For 



both the spin-induced TRS breaking and the increase in 
spin relaxation time we find excellent quantitative agree- 
ment with recent theory. We also show that the quanti- 
tative analysis permits to determine the in-plane g-factor 
of the electrons. 

The three InAlAs/InGaAs/InAlAs quantum wells used 
in our work are very similar to those described in detail 
elsewhere 4.]- Here, we recall that each well is designed 
to have a different (Rashba) SOI strength. The charac- 
teristic spin-split energy A for the different samples is 
A K, 0.5,1.5 and l.SmeV (in what follows we will re- 
fer to these samples as to samples 1, 2, and 3, respec- 
tively). The electron density and mobility at Vgate — OF 
are n ~ 7- lO^^m"^ and /i ~ 4 rn^ /Vs. All measurements 
have been performed on (20 x 80 ^m) Hall-bar shaped 
devices, at 1.6K. A 14 T superconducting magnet is used 
to generate B\\ and home-made split coils mounted on 
the sample holder are used to independently control the 
perpendicular field (_Bx ) • No significant difference in the 
results is observed when the in-plane field is applied par- 
allel or perpendicular to the direction of current flow. 

To understand how an in-plane magnetic field affects 
the electronic transport, we first discuss the behavior of 
sample 1 with the weakest SOI strength. Figure 1 shows 
the magnetoconductance of this sample measured as a 
function of B±^ [2|, for different fixed values of the in- 
plane field B\\. For small values of B\\ (main panel), the 
conductance exhibits a maximum at B±^ = 0, due to 
weak-antilocalization (WAL) superimposed on the back- 
ground of weak- localization (WL) [Sj. As i3|| is increased, 
the amplitude of this maximum is reduced and eventually 
disappears around Bn = 300mT. A further increase in 
Sjl does not result in additional changes of the magneto- 
conductance until B\\ reaches approximately IT (fig. lb). 
Upon increasing B\\ even further, the WL signal is also 
suppressed on the scale of several (~ 10) Tesla (fig.lc). 

These observations allow us to conclude that the sup- 
pression of WAL and of WL in a parallel field are due 
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FIG. 1: The magnetoconductance a{B±) of sample 1 at difFer- 
ent values of By . Three regimes can be identified: increasing 
_B|| from to 350mT results in a suppression of the WAL peak 
(a), increasing _B|| further (up to about By =1T) does not in- 
duce additional changes in the CT(_B_L)-curves (b), for values 
of -B|| larger than IT the WL is suppressed (c). 



FIG. 2: The empty circles are the measured magnetoconduc- 
tance Act = a{B±) — ct(0) of sample 2 at different fixed values 
of _B|| (offset for clarity). By is increased from to IT, in steps 
of O.IT (top to bottom). The solid lines represent best fits to 
the ILP theory. The inset shows the amplitude of the WAL 
peak at B± = as function of _B||, i.e. a{B± — 0, i3||) — ct(0, 0), 
and the best fit to the theory (solid line). 



to two distinct mechanisms causing time reversal syra- 
metry breaking. At large fields, iJy ^ IT, WL (which is 
not sensitive to the spin degree of freedom) is suppressed 
due to TRS breaking caused by the coupling of i3|| to 
the orbital motion of the electrons, owing to the finite 
thickness of the quantum well and the asymmetric con- 
fining potential y. The suppression of the WAL peak 
at smaller values of i?|| originates from a spin-induced 
TRS breaking due to the interplay between By_^(Zeeman 
coupling) and SOI, as predicted theoreticallvjinj . In this 
paper we will focus on the spin mechanism for TRS, and 
discuss the orbital mechanism elsewhere. 

The complete separation of spin and orbital TRS 
breaking, which is essential for the work presented here, 
has not been previously reported Q- In our samples, this 
separation is due to the small QW thickness (« lOnrn) 
and the small effective mass (m* « 0.041mo) which 
make the subband splitting in the QW relatively large, 
as well as to the relatively large gyromagnetic ratio 
[g — 3) [3j ll^ ■ It allows us to account for the magne- 
toconductance curves cr{B±) measured at i3|| < IT in 
terms of existing theories that only consider the coupling 
of i3|| to the electron spin. Therefore, the number of pa- 
rameters that need to be introduced for the quantitative 
analysis of the data is the smallest possible. This makes 
it possible to extract the values of the phase coherence 
time and the spin relaxation time as a function of _B|| with 
great accuracy, as it is needed to observe the dependence 
of Ts on the in-plane magnetic field. 

We have performed a quantitative analysis of the mag- 
netoconductance curves on all samples and for different 
values of n, by fitting the cr{B±) curves with the theory of 



lordanskii, Lyanda-Geller and Pikus (ILP) [Tl[l3- This 
is appropriate for our samples, in which the spin relax- 
ation is governed by the Dyakonov-Percl mechanism [l3j. 
From this analysis, namely from the fits of a{B±) curves 
measured at different values of the in-plane field, we find 
the -By-dependence of Tg and r^, i.e. Ts{B\\) and T0(i3||) 
[lj|. It is worth noting that in the ILP theory only one 
parameter is needed to account for the spin relaxation, 
since Ts{0) = t^JO) = Ts„(0) = 2x^^(0). In the pres- 
ence of an in-plane field, however, these relations may 
not hold, since relaxing the spin along B^\ costs energy 
(« (7/ii3||) whereas relaxation in the direction perpendic- 
ular to B|| does not. Nevertheless, for sufficiently small 
S|j (.g/-ti?|| < kT), the ratios between the different relax- 
ation times are expected to change only minorly under 
the conditions of our experiments. This allows us to treat 
Ts(-B||) as a single fitting parameter. 

Figure 2a displays the results of the fitting procedure 
on sample 2 with the intermediate SOI strength. The 
continuous lines superimposed on the data represent the 
best fit to the ILP theory, and show that the agreement 
between data and theory is excellent for all values of 
Bu . Similar agreement is obtained for the other sam- 
ples and for all the different values of the electron den- 
sity 71. The values of T0(i3||) and rs(i3||), as extracted 
from the fits, are shown in Figs. 3 and 4. Note that, 
since the electron mobility depends on the density, we 
are able to investigate how changing the elastic scatter- 
ing T affects the i3|| -dependence of the phase coherence 
and of the spin-relaxation time. This is of particular 
interest as both r0(_B||) and Ts(-B||) are predicted to de- 







FIG. 3; The symbols represent r^ as a function of -By , as ex- 
tracted from the analysis of the magnetoconductance of sam- 
ple 2, using the ILP theory (see Fig. 2). Different curves cor- 
respond to different values of n (and elastic scattering time 
r). The solid lines are best fits based on the theory describing 
spin-induced dephasing [lO| . The decrease of r^ with decreas- 
ing electron density is consistent with dephasing originating 
from electron-electron interaction. The inset shows the ex- 
tracted T^{B\\) and theoretical fits for sample 1. 
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FIG. 4: The symbols represent Ts as function of By , as ex- 
tracted from the analysis of the magnetoconductance of sam- 
ple 1 (open symbols) and sample 2 (filled symbols). Each set 
of symbols corresponds to a different value of the Dyakonov- 
Perel spin-relaxation time Ts(0), with l/rs(0) — t^r I2T? 
(controlled by changing the gate voltage). The solid lines are 
best fits to the theory (Eq. 3), with k as an added parameter 
(see text). Note that the color and symbol code used in this 
figure for r^ corresponds to that used in Fig. 3 for r^. 



pend on the Dyakonov-Perel spin relaxation time Ts(0) 
(see Eqs. 2 and 3), which is related to r by the relation 
1/t,(0) = t^^TlIT? HI. 

For all values of n, the measured t^{B\\) decreases as 
a function of B^ (Fig. 3), which shows quantitatively 
the breaking of TRS due to the interplay of Zeeman cou- 
pling and SOI. This interplay is predicted to result in a 
quadratic dependence of r^ on B\\ [lO| : 
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where c is a constant given by: 

c = r40)r,(0)(5^MBA)' 



(1) 



(2) 



and gZ is the in-plane g-factor. The solid lines in Fig. 3 
are best fits to the data using Eq. (1) and treating c as a 
(density dependent) fitting parameter. Also in this case 
the agreement between experiment and theory is excel- 
lent for all values of n and for the different samples (the 
inset of Fig. 3 shows the behavior of sample 1. Equally 
good agreement is found for sample 3). 

Using the value of c obtained from fitting the data of 
Fig. 3 we directly obtain gjj (Eq. 2). We find that, 
for each sample, the in-plane g-factor is approximately 
constant as a function of the electron density. The ab- 
solute values are determined to be \gV\ = 2.8 ± 0.1, 
\g*^ = 3.3±0.1and|5||| = 3.5±0.1, for samples 1, 2 and 3, 



respectively. Theoretically, the g-factor in our quantum 
well is predicted to depend substantially on its thickness, 
and is calculated to be \gV\ — 2.8 and \gV\ = 3.5 for 
a thickness of lOnm and 15nm, respectively [l5j- This 
agreement with theory gives additional support of our 
analysis in terms of spin-induced dephasing only, and 
shows that the measurement of WAL in the presence of 
an in-plane field permits to determine the in-plane g- 
factor. Contrary to other methods based on transport 
measurements, this method to determine the g-factor is 
suitable for disordered systems. 

A different way to obtain t^(B\\) (and c), apart from 
fitting the whole cr{B±) curves measured at fixed B\\, 
is by looking at the conductance at B± = as func- 
tion of i?|| . Specifically, the theory for spin-induced de- 
phasing predicts that cr{B± = 0, 0) — a{B± = 0, i?||) = 
filn(T40)/r4B||)) = filn(l + cB2) ^, Also in this 
case, the agreement between theory and data is excellent 
(Fig. 2, inset) and the fitting procedure gives values for 
the parameter c identical to those obtained above. This 
shows the consistency of our quantitative analysis and 
confirms once more the validity of the interpretation of 
the data in terms of spin-induced TRS breaking only. 

Finally, Fig. 4 shows the behavior of the measured 
spin relaxation time as function of Sy for different den- 
sities and different strength of SOI interaction (samples 
1 and 2). In all cases, the measured spin relaxation time 
increases quadratically with increasing the applied in- 



plane field. This directly shows that the presence of an 
in-plane field reduces spin-randomization. The increase 
in rs(i?||) is more pronounced for a small strength of the 
SOI interaction and for short values of the elastic scatter- 
ing time r, i.e. for long Dyakonov-Perel spin-relaxation 
times Ts{0). This is because the Zeeman energy gfiB^^, 
that drives the alignment of the electron spin along Bu , 
competes with the characteristic energy associated to the 
spin-randomization H/ts (0) . 

A quantitative analysis of the data requires a compar- 
ison with theory. For the case of a magnetic field normal 
to the conduction plane, extensive theoretical analysis 
exists [lal- For the case of an in-plane field, however, 
only the relaxation time of the z-component of the spin 
has been calculated as a function of B|| [17]. When the 
Zeeman energy g/i-Bii is much smaller than h/Ts{0), this 
quantity is given by 
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Although theoretical predictions for Ts^(_B||) and 
rsj^(i3||) are not available, we expect Ts^(i?||)/Ts^(0) and 
Ts^ {B^^)/ts^ (0) to exhibit the same functional dependence 
as Ts^{B\\)/ts^{0) as long as (7/X-B|| ^ h/Ts{0) and kT. 
This allows us to compare the measured Ts{B^^)/ts{0) to 
Eq. 3. All the quantities that appear in Eq. 3 are known 
from the previous analysis, and we add a parameter k to 
achieve best fits to the data (theory [l3| predicts k = 1 in 
Eq. 3). Figure 4 shows that in all cases good agreement 
is obtained with k ~ 1 (continuous lines). We conclude 
that the qualitative behavior of the spin-relaxation time 
as a function of i?|| , r and A (or, equivalently, r^ (0)) is the 
one expected, and that, within a small correction factor, 
our results are in quantitative agreement with theoretical 
predictions. 

In view of the quantitative agreement between theory 
and data obtained throughout this work, it is worth con- 
sidering the origin of the small correction factor k. k ^ 1 
may originate from the limited accuracy with which the 
quantities in Eq. 3 are determined. The largest uncer- 
tainty comes from gT, and is approximately 10%. An ad- 
ditional possibility is the i?|| -induced anisotropy of the in- 
plane spin relaxation times, i.e. i?|| breaks spin-rotational 
symmetry in the 2D plane. Although this anisotropy is 
expected to be small for 5J|M-S|| ^ ^/ts{0) and fc^T, as 
mentioned before, it may result in a deviation from k = 1. 
Finally, for sample 1 with the weakest Rashba SOI, the 
Dresselhaus term may not be entirely negligible J18l | . 

In conclusion, we have observed how the partial align- 
ment of the electron spin along an applied in-plane mag- 
netic field determines the orbital and spin dynamics of 
electrons in Rashba 2DEGs. This alignment results in a 
spin-induced time reversal symmetry breaking and in a 
quadratic increase of the spin-relaxation time. The de- 
tailed quantitative analysis of our results demonstrates 



the validity of the existing theory and gives indications 
to its limits. 
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